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Wc present an updated measurement of the CP-odd fraction and the timc-depcndcnt CP asym- 
metries in the decay B° D*+ D*' using (383 ± 4) x 10*^55 pairs collected with the BABm 
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detector. We determine the CP-odd fraction to be 0.143 ± 0.034(stat) ± 0.008(syst). The time- 
dependent CP asymmetry parameters are determined to be C+ = —0.05 ± 0.14(stat) ± 0.02(syst) 
and S+ = —0.72 ± 0.19(stat) ± 0.05(syst). The non-zero value of the measured S+ indicates the 
evidence of CP violation at the 3.7 a confidence level. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



In the Standard Model (SM), CP violation is described 
by a single complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix, V Mea- 
surements of CP asymmetries by the BaBAR Q and 
Belle collaborations have firmly established this ef- 
fect in the b — > (cc)s charmonium decays Q and 
precisely determined the parameter sin2/3, where p is 
a.Tg[-VcdV*^/VtdV^l]. The amplitude of the decay B° 

is dominated by a tree-level, color-allowed h — > 
ccd transition. Within the framework of the SM, the CP 
asymmetry of B^ D*^D*^ is equal to sin2/3 when 
the correction due to penguin diagram contributions is 
neglected. The penguin-induced correction to the CP 
asymmetry, estimated in models based on the factoriza- 
tion approximation and heavy quark symmetry, is pre- 
dicted to be about 2% [5] , while contributions from non- 
SM processes may lead to a large correction |6]. Such 
a deviation in the sin2/3 measurement from that of the 
B° — > {cc)K^*^^ decays would be evidence of physics be- 
yond the SM. 

Studies of CP violation in B" £)(*)±i:)(*)T tran- 
sitions have been carried out by both the BaBAR and 
Belle collaborations. Most recently, the Belle collabo- 
ration reported evidence of large direct CP violation in 
B" ^ D+D- where Cd+d- = -0.91 ± 0.23 ± 0.06 [Z], in 
contradiction to the SM expectation. However, a large 
direct CP violation has not been observed in this chan- 
nel by BaBAR [^], nor in previous measurements with 
B'^ D*'^D*~ decays that involve the same quark-level 
weak decay [l^ . 

The B° D*+D*- decay proceeds through the CP- 
even S and D waves and through the CP-odd P wave. 
In this Letter, we present an improved measurement of 
the CP-odd fraction R± based on a time-integrated one- 
dimensional angular analysis. We also present an im- 
proved measurement of the time-dependent CP asymme- 
try, obtained from a combined analysis of time-dependent 
flavor-tagged decays and the one-dimensional angular 
distribution of the decay products. 

The data used in this analysis comprise (383 ±4) x 10^ 
T{AS)-^ BB decays collected with the BaBAR detec- 
tor [ll| at the PEP-II asymmetric-energy e+e^ storage 
rings. We use a Monte Carlo (MC) simulation based on 
GEANT4 [13] to validate the analysis procedure and to 
study the relevant backgrounds. 

We select B^ — > D*^D*~ candidates from oppositely 
charged pairs of D* mesons. The is reconstructed 
in its decays to D'^tt'^ and ZJ+tt". We reconstruct candi- 
dates for and £)+ mesons in the modes Z?" K^tt^, 



K'TT+TT^, K-tt+tt+t:-, K'^tt+tt- and D+ K-tt+tt+. 
We reject the _B" candidates for which both D* mesons 
decay to Dtt'^ because of the smaller branching frac- 
tion and larger backgrounds. To suppress the e^e^ — > 
qq [q = u, d, s, and c) continuum background, we require 
the ratio of the second and zcroth order Fox- Wolfram 
moments [l3| to be less than 0.6. 

For each B'^ D*^D*^ candidate, we construct a 
likelihood function £mass from the masses and mass un- 
certainties of the D and D* candidates |1^]. In this like- 
lihood, the D mass resolution is modeled by a Gaus- 
sian function whose variance is determined candidate- 
by-candidate from the mass uncertainty resulting from a 
vertex fit of the D meson decay products. The D* — D 
mass difference resolution is modeled by the sum of two 
Gaussian distributions whose parameters are determined 
from simulated events. The maximum allowed values of 
— ln£mass and \AE\ = — iSbeamI, the difference be- 
tween the 5° candidate energy and the beam energy 
Ehcam in the T(4S') frame, are optimized separately for 
each final state using simulated events to obtain the high- 
est expected signal significance. 

We include candidates with an energy-substituted 
mass, = V ^bcam ~ P*B ' greater than 5.23 GeV/c^, 
where p*^ is the candidate momentum in the T(45') 
frame. On average, we have 1.8 B° candidates per event 
in data after all the selection requirements. In cases 
where more than one candidate is reconstructed in an 
event, the candidate with the smallest value of — In/Imass 
is chosen. Studies using MC samples show that this pro- 
cedure results in the selection of the correct 5° candidate 
more than 95 % of the time. 

The total probability density function (PDF) of the 
Toes distribution is the sum of the signal and background 
components. The signal PDF is modeled by a Gaussian 
function and the combinatorial background is described 
by a threshold function Studies based on MC sim- 
ulation show that there is a small peaking background 
from B^ — > D*'^D*^ in which a originating from a 
_D*° decay is combined with a random soft tt" to form 
a D*~ candidate. This background is described by the 
same PDF as the signal, and its fraction with respect to 
the signal yield is fixed to (1.8 ± 1.8)%, as determined in 
MC simulation. An unbinned maximum likelihood (ML) 
fit to the 771ES distribution yields 617 ± 33(stat) signal 
events, where the mean and width of the signal Gaussian 
function and the threshold function shape parameters are 
allowed to vary in the fit. The signal purity in the region 
of Toes > 5.27 GeV/c^ is approximately 65%. 
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Following [1^, we define three angles depicted in Fig- 
ure [1] within the transversity framework: the angle 9i 
between the momentum of the slow pion from the D*^ 
and the direction opposite the D*~^ flight in the D*~ rest 
frame; the polar angle 0tr and azimuthal angle (jjtr of the 
slow pion from the D*'^ evaluated in the D*+ rest frame, 
where the coordinate system is defined with the z axis 
normal to the D*~ decay plane and the x axis opposite 
to the D*~ momentum. 




FIG. 1: Depiction of the B° D*+ D*' decay in the 
transversity basis. The three transversity angles are defined 
in the text. 

The time-dependent an gula r distribution of the decay 
products is given in Ref. [l7|. Taking into account the 
detector efficiency as a function of the transversity angles 
and integrating over the decay time and the angles 9i and 
4>f[, we obtain a one-dimensional differential decay rate: 



1 dV 



r d cos 0t 
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1 + a 



{(l-i?±)sin2 0tr 

/o(cOS ^tr 



1 



where R±^ = |Aj 
l^||n/(l^ol^ + 



2 ' 2 

2R±_ cos^ 9t^ X /^(cos6'tr)} , 

|V(|Ao|2 + |A|||2 + |Axn,a 



/||(cos 6't,) 



(1) 



(l^ol 



Aqi is the amplitude for longi- 



tudinally polarized D* mesons, A\\ and A±^ are the 
amplitudes for parallel and perpendicular transversely 
polarized D* mesons. The three efficiency moments 
/jt(cos ^tr), where (fc — 0, ||,-L), are defined as 

Jfe (cos 6*11) = J dcoseid(l)tr 9ki6i,M£i6i,dti-,M, (2) 

where 50 = 4 cos^ 9i cos^ 0tr, 5|t = 2 sin^ 9i sin^ (j)tr, gj_ = 
sin^ 9i , and e is the overall detector efficiency. The effi- 
ciency moments are parameterized as second-order even 
polynomials of cos 0tr with parameter values determined 
from the MC simulation. In fact, the three Ik functions 
deviate only slightly from a constant, making the decay 
distribution (Eq. [T|) nearly independent of the amplitude 
asymmetry a. 



The CP-odd fraction R± is measured in a simulta- 
neous unbinned ML fit to the cos^tr and the tties dis- 
tributions shown in Figure [51 The background in the 
cos 9tr distribution is modeled as an even, second-order 
polynomial, while the signal PDF is given by Eq. [TJ 
The finite detector resolution of the 9ti- measurement is 
modeled by the sum of three Gaussian functions plus a 
small tail component that accounts for misreconstructed 
events, where all the parameters are fixed to the val- 
ues determined in the MC simulation. The resolution 
function is convolved with the signal PDF in the max- 
imum likelihood fit. We categorize events into three 
types: D^+D*- {D^tt+ ,Wtt-), {D°Tr+ , D-7t°), and 
(ZJ+TT^, ZJ^TT"), each with different signal-fraction param- 
eters in the likelihood fit. Their efficiency moments and 
cos ^tr resolutions are separately determined from the MC 
simulation. The other parameters, determined in the 
likehhood fit, are the cos^tr background-shape parame- 
ter, three toes parameters (width and mean of the signal 
Gaussian, and the threshold function shape parameter), 
as well as i? I . 
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FIG. 2: Measured distribution of mES (a) and of cos^tr in the 
region tties > 5.27 GeV/c^ (b). The solid line is the projection 
of the fit result. The dotted line represents the background 
component. 

After fitting to data and taking into account possible 
systematic uncertainties, we find 



Ri_ = 0.143 ± 0.034(stat) ± 0.008(syst) . 



(3) 



Figure [2] shows the projections of the data and the fit 
result onto toes and cos6'tr- 

In the fit described above, the value of a, the asymme- 
try between the two CP-even amplitudes in the transver- 
sity framework, is fixed to zero. We estimate the cor- 
responding systematic uncertainty by varying its value 
from —1 to +1 and find negligible change (0.003) in the 
fitted value of R±. Other systematic uncertainties arise 
from varying fixed parameters within their errors: the 
parameterization of the angular resolution (0.006), the 
determination of the efficiency moments (0.004), and the 
background parameterization (0.004). The total system- 
atic uncertainty on R± is 0.008. 
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We perform a combined analysis of the cos^tr dis- 
tribution and its time dependence to extract the time- 
dependent CP asymmetry, using the event sample de- 
scribed previously. We use information from the other 
B meson (i3tag) in the event to tag the initial flavor of 
the fully reconstructed —> D*+D*^ candidate (Brcc)- 
The multivariate flavor tagging algorithm is described in 
detail elsewhere [3]. We define six mutually exclusive 
tagging categories in order of expected tag purity from 
lepton to hadron, which includes kaon and pion tags. 
The total effective tagging efficiency of this algorithm is 
(30.5 ± 0.4) %. 

The decay rate /+(/_) for a neutral B meson accom- 
panied by a B^{B^) tag is given by 

/±(At,cos6't,.) oc e^l'^*l/^««|G(l T Aw) ± (1 - 2uj) 

[Fsin(Am<jAt) - cos(ATOdAi)]|, (4) 

where Ai = t^-^c — ^tag is the difference between the 
proper decay time of the i^rcc and i?tag mesons, r^o = 
(1.530 ±0.009) ps is the B^ lifetime, and Am^ = (0.507± 
0.005) ps~^ is the mass difference between the B^-B^ 
mass eigenstates [l^. The average mistag probability u) 
describes the effect of incorrect tags, and Au; is the dif- 
ference between the mistag rate for and B^ . The G, 
F and H coefficients are defined as: 

G = (l - R.i_)sm^ 9tr + 2R^cos'^ 9tr, 

F = {l-R.j_)S+sm'^9tr-2R^S^cos'^0tr, (5) 

H = {l-Rj_)C+ sin^ 0t, + 2R^C± cos^ Ot,, 

where we allow the three transversity amplitudes to have 
different Xk = {q/p){Ak/Ak) {k = 0, ||, 1) [17] due to pos- 
sibly different penguin-to-tree amplitude ratios, and de- 
fine the CP asymmetry parameters Gfc = (1 — |Afcp)/(l-|- 
|Afep), Sk = 2Jm{Xk) / {I + |Afep). Here, we also define: 

^ qil^llP + Col^op ^ S|||A||p + 5o|Aop 
^+ |^||P + |AoP ' + |A||P + l^oP • ^ ' 

In the absence of penguin contributions, we expect that 
Co = q, =C±^0, and 5*0 = 5|| = 5_l = - sin2/3 ^. 

In Eq. m the small detector efficiency effects are not 
taken into account and instead are absorbed into the 
value of R±^ which is allowed to vary in the final fit. 
Any bias in the resulting values of C+, C_l, 5*+, and S± 
is below the sensitivity of our MC validation sample and 
is accounted for in the MC statistics systematic. Hence, a 
dedicated method to correct for detector efficiency is not 
required. However, the "effective" value of R± obtained 
in this way is not identical to the value measured from 
the time-integrated analysis that includes the efficiency 
correction. This approach incorporates the uncertainty 
in R± directly into the determination of the CP param- 
eters in the ML fit. 



The technique used to measure the CP asymmetry is 
analogous to that used in BaBAR measurements as de- 
scribed in Ref. We calculate the time interval 
At between the two B decays from the measured sep- 
aration Az between the decay vertices of Brec and i^tag 
along the colHsion (z) axis [2^. The z position of the B^ec 
vertex is determined from the charged daughter tracks. 
The i3tag decay vertex is determined by fitting charged 
tracks not belonging to the Brec candidate to a common 
vertex, employing constraints from the beam spot loca- 
tion and the i?rec momentum (20| . Only events with a 
At uncertainty less than 2.5 ps and a measured \At\ less 
than 20 ps are accepted. We perform a simultaneous un- 
binned ML fit to the cos^tr, At, and rngg distributions 
to extract the CP asymmetry. The signal PDF in 6'tr and 
At is given by Eq. [D The signal mistag probability and 
the difference between the mistag rate for 5° and are 
determined for each tagging category from a large sam- 
ple of neutral B decays to flavor eigenstates, i?fiav. In the 
likelihood fit, the expression in Eq. [4]is convolved with 
an empirical At resolution function determined from the 
-Bflav sample. The 6tr resolution is accounted for in the 
same way as described previously. 

Our increased statistics allows for better treatment of 
the background in this analysis. The background At dis- 
tributions are parameterized by an empirical description 
that includes components that have zero lifetime, and 
that have an effective lifetime similar to the signal. The 
lifetime of the second component and its relative fraction 
are allowed to vary in the likelihood fit. We also allow the 
lifetime component to have free effective CP asymmetry 
parameters, Ccff and S'efi, for each tagging category to 
take into account a possible difference in mistag rates in 
the background. The background shape in 6'tr is modeled 
by an even, second-order polynomial in cos On , as in the 
time- integrated angular analysis. 

From our fit to data we determine 

C+ = -0.05±0.14(stat)±0.02(syst), 

C± = 0.23±0.67(stat)±0.10(syst), 

S+ = -0.72±0.19(stat)±0.05(syst), 

Si_ = -1.83±1.04(stat)±0.23(syst). (7) 

The correlations between C+ and C_ and between S+ 
and S- are —0.46 and 0.39 respectively. All other corre- 
lations are negligible. Figure[3]shows the At distributions 
and asymmetry in yield between B'^ and B^ tags, overlaid 
with the result of the likelihood fit. Because R± is small, 
we have rather large statistical uncertainties for the mea- 
sured C± and S± values. We repeat the fit assuming 
that both CP-even and CP-odd states have the same CP 
asymmetry, i.e. C+ = C± = C and S-^ = S± ^ S. We 
find 

C -0.02±0.11(stat) ±0.02(syst), 

S = -0.66±0.19(stat)±0.04(syst). (8) 
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In both cases, the effective CP asymmetries in the back- 
ground are found to be consistent with zero. To further 
test the consistency of the fitting procedure, the same 
analysis is apphed to the i?° D*~^D*~ control sample. 
The result is consistent with no CP violation as expected. 




4 6 8 
At(ps) 

FIG. 3: The distribution in At of the yield in the region 
rriES > 5.27 GeV/c^ for B° tagged candidates (a) and 

the raw asymmetry {Ngo — Ngo)/ (Ngo + N^a), as functions 
of At (b). In (a), the solid (dashed) curves represent the fit 
to the data for B° (5°) tags. 



The sources of systematic uncertainty on the CP asym- 
metries and their estimated magnitudes are summarized 
in Table HI We vary the yield and CP asymmetries of 
possible backgrounds that peak under the signal. We 
also vary fixed parameters in the fit for the assumed pa- 
rameterization of the At resolution function, the possible 
differences between the -Bflav and Bcp mistag fractions, 
and knowledge of the event-by-event beam-spot position. 
We evaluate the uncertainty due to possible interference 
between the suppressed b — > ucd amplitude and the fa- 
vored b 



cud amplitude for some tag-side decays 21 1. 
We also include systematic uncertainties incurred from 
the finite MC sample used to verify the fitting method. 
All of the systematic uncertainties are much smaller than 
the statistical uncertainties. 

In summary, we have reported measurements of the 
CP-odd fraction, _Rx, and time-dependent CP asymme- 
tries for the decay i?° D*~^D*~. The measurement 
is consistent with and supersedes the previous BaBAR re- 
sult The time-dependent asymmetries are found to 
be consistent with the SM predictions. The non-zero 
value of the measured S+ indicates the evidence of CP 
violation at the 3.7 cr confidence level. 
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Source 


c+ 


s+ 




Si_ 


C 


s 


Peaking backgrounds 


0.008 


0.028 


0.037 


0.110 


0.003 


0.028 


At resolution parameterization 


0.009 


0.011 


0.018 


0.022 


0.008 


0.010 


Mistag fraction difTerences 


0.008 


0.024 


0.016 


0.035 


0.008 


0.024 


Beam-spot position 


0.004 


0.007 


0.019 


0.042 


0.003 


0.005 


Arud, TB 


n nn a 
0.004 


O.OOb 


O.OlD 


0.004 


0.001 


O.OOb 


Angular resolution 


0.009 


0.031 


0.076 


0.116 


0.008 


0.012 


Tag-side interference and others 


0.014 


0.009 


0.017 


0.021 


0.014 


0.009 


MC statistics 


0.005 


0.013 


0.031 


0.150 


0.001 


0.013 


Total 


0.024 


0.053 


0.098 


0.229 


0.021 


0.044 



TABLE I: Systematic errors on time-dependent CP asymmetry parameters for the decay B° — » . 



8 



